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Abstract: We briefly present the super symmetric double Darboux method and next apply it 
to the continuum of the quantum Taub cosmological model as a toy model in order to generate 
a one-parameter family of bosonic Taub potentials and the corresponding wavefunctions. 

Recently, extensive work has been devoted to generating quasi-isospectral potentials by the 
fj^ ■ methods of supersymmetric quantum mechanics (SUSYQM) [|J . This work starts from the one 

dimensional (ID) Schrodinger equation with a given potential (called "bosonic" potential V_), 
whose "ground-state" (nodeless) wave function is known too. The method allows one to generate 
families of new potentials, which may look quite different form the original one, but have the 
. same discrete spectrum up to the ground state. SUSYQM is in fact based on a procedure already 

CN | known to Darboux We shall use an extension of SUSYQM formalism based on the general 

Riccati equation and on a state Uq(x) of a potential V_(x), by which one can generate a one- 
. parameter family of bosonic potentials V(x; A) including V_(x) |J. In SUSYQM, the labeling, 

real parameter A should be positive for radial problems and should not be in the interval [—1,0] 
for the full line problems, but in quantum cosmology this constraints can be relaxed. 
We start with the ID Hamiltonian 

a* d 2 

H = "^2+V-(X), (1) 



which can be factored as H = A + A by means of the operators A and A + defined as 



A=A + W ( X ) a+ = ~ + W(x). (2) 

dx dx 

In this way, one can get V_(x) = W 2 (x) — , where W(x) is called the superpotential. If 
Uq{x) is the ground state solution of the ID Schrodinger equation, the superpotential is given 
by 

W(x) = -A[ lnZ4(x )]. ( 3 ) 

In the supersymmetric scheme there exists another potential, called "fermionic" potential, which 
is the partner of the bosonic one, that can be written as 



dx dx 2 



V+(x) = W 2 + — =V_-2- T [lnW (x)]. (4) 



It corresponds to the supersymmetric partner Hamiltonian H' = AA + = — 4-? + V+ (x) . 
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If the potential V_(x) has eigenfunctions U n {x) at energies E n , then the SUSY partner 
potential V+(x) has the same energy eigenvalues as V_(x) with eigenfunctions AW n (x), except 
that there is no ground state at E=0 since AUo(x) = 0. This is the standard procedure for 
deleting the ground state and obtaining V+(x). Now, suppose one asks for the most general 
superpotential such that V+(x) = W 2 + It is easy to see that one particular solution 

to this equation is W p = W(x), where W(x) is given by Eq.@. One is led to consider the 

following Riccati equation W 2 + ^ = W 2 + ^j 2 -, whose general solution can be written down 
as W(x) = W p (x) + where v(x) is an unknown function. Using this ansatz, one obtains for 
the function v(x) the following Bernoulli equation 

dv(X) 2v(x)W p (x) = l, (5) 



dx 



that has the following solution 



Zq(x) + A 
W 2 (x) 

where 1q(x) = J x U^iy) dy. Thus, W(x) can be written as follows 

W(x) = Wp (x) + ^-[ln(l (x) + X)} 



V ( X ) = 7/2/.^ » ( 6 ) 



d ZYp(x) 

"dx [ln( 2o(x)TA )] - (7) 



Finally, one gets easily the family of "bosonic" potentials 

dW(x; A) 



V(x;A) = W 2 (x;A) 



dx 

d 2 



= V_(x)-2 5 - f [ln(J (x) + A)] 

4^ (xK(x) 2^ 4 (x) 
" V - (X) " T (x) + A + (T (x) + A) 2 - (8) 

All V(x; A) have the same supersymmetric partner potential V+(x) obtained by deleting the 
ground state. 

On the other hand, by inspection of Eqs. @ and (|7|) we can obtain the new family wave 
function for the potentials V(x; A) 

In SUSYQM, some authors have argued that in order to have normalizable wave functions, 
the function f(A) must be of the following form f(A) = y A(A + 1). A connection with other 
isospectral methods has been found, by noticing that the limiting values -1 and for the param- 
eter A lead to the Abraham-Moses procedure [|J| and Pursey's ||] one, respectively. On the other 
hand, in quantum cosmology we shall work with /(A) = 1, i.e., with nonnormalizable solutions. 

The supersymmetric double Darboux method can be applied to any one-dimensional system, 
whose dynamics is dictated by a Schrodinger-like equation. Here we shall apply it by employing 
a wave function belonging to the continuum of the quantum minisuperspace Taub model, which 
is one of the simplest models in quantum cosmology. It has been studied in some detail by 
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Ryan and collaborators ||], who found that it is separable. Indeed, the Taub Wheeler-DeWitt 
equation is 

d 2 ^ d 2 ^ 

___ + e -V(ffl* = 0, (10) 

where V(/3) = ^(e _8/3 — 4e -2 ^). Eq.(|To|) can be separated in the variables x\ = Aa — 8/3 and 
X2 = 4a — 2/3. Thus, one gets two independent ID problems for which the supersymmetric 
procedures are standard practice. The two equations are as follows 

, 1 X! u 2 

~^ +1 r X2u2 = wV (12) 

The quantity co is mathematically the separation constant, which physically is related to the 
wavenumber of a positive energy level. 

Martinez and Ryan have considered a wavepacket solution made of wavefunctions having 
the form of a product of modified Bessel functions of imaginary order. We shall slightly modify 
their ^ as follows 

* = Ul u 2 = KU\^ 1/2 )[UU\^ 2 ' 2 ) + W^ 2 )] (13) 
b 3 3 

since, according to Dunster []?]], the L function defined as 

7ri 

L 2 iu, = n . , /r< r{hiu) + I-2iw) (14) 

2smh(2u;7r) 

contrary to the 1^ function used by Martinez and Ryan, being real on the real axis is a better 
companion for the K function of imaginary order, which is known to be real on the real axis. 

In the following, we shall make the double Darboux construction on the base of a ^ wave- 
function of the type given in Eq. (|l3|). 

The family of bosonic potentials for the x\ variable will be 



V^Ai) = ^e xi - 2^2 ln[Ai + J( Xl )] (15) 



and for the x 2 one 
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V 2 (x 2 ; A 2 ) = -e X2 - 2—* ln[A 2 + J(x 2 )] (16) 



dx^ 

where the integrals are 

/ x i 1 /' x 2 2 2 

KL(^/ 2 )dy, X(x 2 ) = / [W-e^ 2 ) + K 2iw (-e^ 2 )] 2 dy . (17) 

The total wavefunction of the bosonic family of Taub potentials has the following form 
AT, x K^(ie^) (L 2i ^(|e^) + K 2iL ,(|e^)) 

* (xi ' X2;Ai ' A2)=UiU2= [x(x i)+ a 1 ] m^Tx;] • (18) 

In conclusion, we used the SUSYQM double Darboux method to introduce a one-parameter 
family of "quantum" Taub potentials having a different set of wavefunctions in the continuum 
region of the spectrum with respect to the original potential. For more details see ref. B. 
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The double Darboux method appears to be a quite general and useful method to generate 
new sets of quantum cosmological solutions for the one-dimensional cases. This is so because 
any potential in the Schrodinger equation has a classical continuum of positive energy nonnor- 
malizable solutions. Ryan and collaborators || were among the first to pay attention to the 
continuum part of the Wheeler-DeWitt spectrum. Selecting by means of some preliminary phys- 
ical arguments that are of quantum scattering type one of the continuum solutions ||], one can 
perform the double Darboux construction on that solution and generate by this means families of 
cosmological potentials as well as strictly isospectral cosmological wavefunctions. The quantum 
Taub model just illustrates this nice feature of the double Darboux method. The parameter A 
looks like a decoherence parameter embodying a sort of cosmological dissipation (or damping) 
distance. Its true physical origin, i.e., the right interpretation, should be further investigated. 

Finally, one can easily apply combinations of any pairs of Abraham-Moses procedure, Pursey's 
one, and the Darboux one. However, only the double Darboux method used here leads to re- 
flection and transmission amplitudes identical to those of the original potential. 
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